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Abstract 

We employ an isospin dependent version of the QMD transport model to study the influence of the isospin 
dependent part of the nuclear matter equation of state and in-medium nucleon-nuelcon cross-sections on the 
dynamics of heavy-ion collisions at intermediate energies. We find that the extraction of useful information 
on the isospin-dependent part of the equation of state of nuclear matter from proton or neutron elliptic 
flows is obstructed by their sensitivity to model parameters and in-medium values of nucleon-nucleon cross- 
sections. Opposite to that, neutron-proton elliptic flow difference shows little dependence on those variables 
while its dependence on the isospin asymmetric EoS is enhanced, making it more suitable for a model 
independent constraining of the high-density behaviour of asy-EoS. Comparison with existing experimental 
FOPI-LAND neutron- hydrogen data can be used to set an upper limit to the softness of asy-EoS. Successful 
constraining of the asy-EoS via neutron-proton elliptic flow difference will require experimental data of 
higher accuracy than presently available. 

Keywords: symmetry energy, equation of state of nuclear matter, quantum molecular dynamics, 
heavy-ion collisions, elliptic flow 



1. Introduction 

One of the remaining opened questions in nuclear physics is the equation of state (EoS) of isospin 
asymmetric nuclear matter, i.e. the density dependence of the symmetry energy (SE). Its precise knowledge 
is mandatory for a proper understanding of nuclear structure of rare isotopes, dynamics and products 
of heavy-ion collisions, and most importantly for astrophysical processes such as neutron star cooling and 
supernovae explosions [1, 2]. Intermediate energy nuclear reactions involving radioactive beams have allowed, 
by studying the thickness of neutron skins, deformation, binding energies and isospin diffusion to constrain 
the density dependence of SE at densities below saturation [3, 4]. Existing theoretical models describing 
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its density dependence generally agree with each other in this density regime, but their predictions start to 
diverge well before regions with densities p > 2po are reached [5]. 

Heavy-ion collisions are the only available opportunity to recreate in laboratory fractions of nuclear 
matter with densities significantly above saturation density, presenting the opportunity for a quantitative 
study of the high-density behavior of the symmetry energy together with a few other closely related topics: in- 
medium nucleon-nucleon (NN) cross-sections, in-medium nucleon masses and the EoS of symmetric nuclear 
matter. These topics are aspects of the same, more fundamental, problem: the in-medium NN interaction. 

A quantitative study of the density dependence of SE demands the identification of observables which 
exhibit a large sensitivity to changes in its magnitude, but show small or no change with respect to variations 
of other variables or model parameters (compressibility modulus, elastic nucleon-nucleon cross-sections, etc.). 
Such observables have been a topic of intensive search during the last decade and a couple of promising 
candidates have been identified: the ratio of neutron/proton ratio of squeezed out nucleons [6], light cluster 
emission [7], 7r^/7r+ multiplicity ratio in central collisions [8, 9], double neutron to proton ratios of nucleon 
emission from isospin-asymmctric but mass-symmetric reactions [10], N/Z dependence of balance energy [11] 
and others. The comparison of theoretical predictions for the tt^ /n^ multiplicity ratio with the experimental 
values of the FOPI collaboration [12] has yielded contradictory results: the study in Ref. [8] points towards 
a soft asy-EoS at supranormal densities, while the results of Ref. [9] favor a stiff scenario. Identifying 
observables that can be compared to currently existing or in the near future available experimental data is 
therefore of great importance. 

This work is devoted to the study of other two observables that have been previously shown to present 
an important sensitivity to changes in the isospin dependent part of EoS: elliptic flow (EF) [13, 14, 15] and 
neutron-proton elliptic flow difference (EFD) [16, 17, 18, 19]. Their dependence on the EoS of both symmetric 
and asymmetric nuclear matter, microscopic nucleon-nucleon cross-sections and model parameters is studied 
in detail and exemplified for the case of Au+Au collisions at incident energy of 400 AMeV. A comparison of 
our proton EF estimates with the high precision data of the FOPI experiment is performed [20] . Additionally, 
a comparison of the theoretical values of neutron-proton EFD with the FOPI-LAND experimental neutron- 
hydrogen and neutron-proton EFD data [21, 22] is presented. 

The Letter is structured as follows: in Section II we review the basic ingredients of the QMD transport 
model we have used and then present the pertinent details regarding the in-medium NN microscopical 
cross-section and EoS of nuclear matter we have employed. We continue with Section III in which results 
on proton and neutron EF and neutron-proton EFD of simulated Au+Au collisions are presented with an 
emphasis on their sensitivity to various ingredients of the transport model together with a comparison of 
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available FOPI and FOPI-LAND experimental data. We end with a Section dedicated to conclusions. 



2. The transport model 

2.1. The QMD transport model 

The heavy-ion collision dynamics is described within the framework of quantum molecular dynamics 
(QMD) [23; 24], a semiclassical transport model which accounts for relevant quantum aspects like stochastic 
scattering and Pauli blocking of nucleons. The Tiibingen transport model code [25], which constitutes 
the backbone of the code employed to produce the results reported in this work, has been expanded by 
adding explicit density dependence of the microscopical nucleon-nucleon cross-sections and building in a 
isospin dependent EoS. Two model parameters, the spread of the single nucleon wave function L and the 
compressibility modulus K, are of relevance for the present study. There values have been set to 2L^^8 
fm^ and K~210 MeV if not otherwise stated. Further details of the transport model can be found in 
Refs. [25, 26]. 

2.2. The in-medium nucleon-nucleon interaction 

One of the important ingredients of the transport model utilized for the simulation of heavy-ion collisions 
are the microscopic elastic and inelastic cross-sections. For the case of elastic collisions there are both 
theoretically [27, 28] and experimentally driven [29, 30] studies hint for important in-medium modifications. 
For the case of total cross-sections the following density dependent parametrizations are being used ^ [27, 28] 

a^AEia,,p) = [31.5 + 0.092 X (20.2 -EU'r]xh2±^^^0^ (1) 

a,,iE,^,, p) = [23.5 + 0.256 x (18.501 - Elff-^] x g^pf ^r2^' 

for energies below the pion production threshold; here Eiab is the incident kinetic energy in MeV, while the 
density p is expressed in fm~'^. In-medium modifications of the angular dependence of the elementary np 
and pp differential cross-sections [27, 28] are neglected in the present study as they were shown to have small 
impact on observables [31]. Values of the elementary cross-sections of reactions involving at least one excited 
baryon (either A or N*) arc supposed to remain unmodified in a dense nuclear medium. Alternatively the 
Cugnon parametrization of vacuum NN cross-sections [32] can be used. The only sizable difference between 



^The parametrization of a^p we have used differs slightly from the one in [28] in order to reproduce accurately the theoretical 
estimates for in-medium a^p listed in Table I of that reference; the parametrization proposed in Eq. (1) of Ref. [28] is rather 
inaccurate in that respect. 
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the two mentioned paramctrizations occurs in the neutron-proton channel at incident energies below 100 
MeV. 

The in-medium NN cross-sections should also exhibit a dependence on the isospin asymmetry factor /3. 
Most of existing works on the in-medium NN interactions have been focused on isospin symmetric nuclear 
matter. The authors of Ref. [33] have therefore extended the effective mass scaling model for the in-medium 
NN cross-sections to isospin asymmetric nuclear matter. The central assumption of this model is that the NN 
matrix elements retain their vacuum analytical expressions in the medium, leading to a isospin asymmetry 
dependence indirectly through the expressions of the effective nucleon masses as functions of /? 

m^vi (p, = 0) niN^ [p^p = 0) 

The above expression is used with minimal changes above the pion production threshold to describe both 
the density and isospin asymmetry dependence of NN cross-sections by considering a dependence of nucleon 
mass on both parameters. The in- medium nucleon mass is derived from the momentum dependent part of 
the mean field potential presented in the next section by making use of the expression: mAr(p, /3)/tojv(p — 
0,/3 = 0) = + mN{p = 0,^ = 0)/p X dU/dp). 

2.3. Equation of state of isospin asymmetric nuclear matter 

The central ingredient for the study presented in this Letter is the isospin dependent part of the EoS 
of nuclear matter. Theoretical studies have shown that the EoS of asymmetric nuclear matter can be 
approximately expressed as £(p, /3) = £{p) + £sym{p) 0^ where (3 = {pn — Pp)/{pn + Pp) is the isospin 
asymmetry. The value of the symmetry energy (SE), £syrm is approximately known only at saturation 
density: £sym{po~'27-S6 MeV). For its density dependence, microscopical models lead to divergent results 
at supranormal densities. We will be using a parametrization that has been developed in Ref. [34] starting 
from the Gogny effective interaction in order to obtain an explicit momentum dependence of the symmetry 
energy part. It reads 

[/(p,/3,p,T,x) = A^{x)^+Mx)^+B{p/por{l-x/3')-8rx^^(3pr> (3) 

PO PO Cr + 1 pfj 

2Ctt f j3 , fri^iP') 2CrT' f ,3 , /r'(?', P') 



pa J l + {p-^Y/K^ PO J l + ip-p'r/A^ 

where the parameters Au{x)., Ai{x), B,Ct,t' ,Ct,t and A have been obtained by fitting the momentum- 
dependence of U{p, P,p,T,x) to the one predicted by the Gogny potential and requiring that the properties 
of nuclear matter at saturation density are reproduced. The compressibility modulus K is set to be a soft 
one (210 MeV) and the value of energy symmetry at po is set to 30 MeV. The parameter x can be adjusted to 
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Figure 1: Left Panel: Dependence of the proton elliptic flow parameter ™ AuAu collisions at 400 AMeV on in-medium 
NN cross-sections. The isospin independent EoS is used. Right Panel: Dependence of on the asy-EoS; the vacuum Li- 
Machleidt A'^A'^ cross-sections have been used. Experimental data points are due to the FOPI collaboration [20] . The notation 
"Li-Machlcidt abed" stands for inclusion (1) or omission(O) of nuclcon-nuclcon cross-sections dependence on the following 
parameters: a - density below pion production threshold (ttPT) (paramctrization in Equation (1)), b - isospin asymmetry 
below ttPT, c - density above ttPT and d- isospin asymmetry above ttPT. In the scenarios fe, c and d the scaling law of 
Equation (2) has been applied. 

mimic the density dependence of the symmetry energy of various microscopical theoretical models and has 
been varied in the current study between -2 and 2. The value x = —2 corresponds to a super-stiff asy-EoS, 
X = 2 to a super-soft asy-EoS, while setting x = the momentum independent part of SE is omitted. Most 
theoretical predictions for the density dependence of SE lie in the region delimited by the x~-2 and x—1 
curves [5]. Further details, including the values of the parameters used in the calculations can be found 
in [33, 34]. 

3. Experimental observables 

3.1. Elliptic Flow 

The azimuthal distribution of protons (or neutrons) resulted in heavy-ion collisions can be approxi- 
mately described by dN/d(f) = {N/{2'k) [l -I- vicoscj) + U2Cos2(/)], vi and V2 being called the sidewards 
and elliptic flow (EF) parameters respectively. The elliptic flow can by extracted from simulated or ex- 
perimental data by computing the following average over the respective particle specie in the final state 
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Figure 2: Elliptic flow of neutrons (left panel) and protons (left panel) in Au-f Au eoUisions at 400 AMeV incident energy 
and impact parameter b=5.5 fm for given choices of the NN microscopical cross-sections as a function of the variable x that 
parametrizes the isospin dependence of the asy-EoS. 



We start by comparing, in Figure 1, the simulated values for the proton EF parameter v\ to the ex- 
perimental values obtained by the FOPI collaboration for AuAu collision at 400 AMeV incident energy. 
The results where obtained by imposing the following cuts on the output particle spectra: 80° < (f> < 100° 
and 0.8 < pt/pcm 1^ 1-8; additionally, the system of reference is rotated around the y-axis with an an- 
gle 8fioio [35] along the direction of flow. In the left panel the dependence of the simulated values on 
the NN cross-section is displayed, with the isospin dependence omitted. There is little difference between 
the two parametrizations of the vacuum NN cross-sections, namely the ones due to Cugnon [32] and Li- 
Machleidt [27, 28] respectively. The sensitivity of to the in-medium NN cross-sections is moderate in size 
with the exception of the case when a density dependence above ttPT is considered. This is partially due to 
the fact that the flow angle depends strongly on the scenario used for the in-medium NN cross-section 

and as a consequence much of the EF dependence on in-mediumA'^A'^ cross-sections is rotated away. 

In the right panel of Figure 1 the sensitivity of f 2 to asy-EoS is studied. A monotonous dependence on 
the X parameters is observed with the exception of peripheral collisions in conjunction with a super-soft asy- 
EoS. A similar conclusion holds for neutrons with the exception of an opposite direction of the monotonous 
variation of the EF parameter with increasing x. This behavior is more clearly visible in Fig. 2where the 
values of neutron (left panel) and proton (right panel) EF as a function of the parameter x are plotted for 
mid-central collisions. We notice, additionally, that the difference (a; = 2) — V2{x = —2) is in magnitude 
larger than the observed sensitivity on the NN cross-sections in the left panel of Figure 1. 

The most reliable extraction of the compressibility modulus K of nuclear matter has been made possible 
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by studying the multiplicity ratio of production in heavy (Au+Au) over light (C+C) nuclei at incident 
energies close to 1 AGeV by the KaoS Collaboration [36, 37, 38] pointing towards a soft EoS. At lower 
incident energies the situation is not as clear: KaoS result points to an even softer EoS while the study of 
sidewards flow by the FOPI collaboration [39] favors a soft or hard EoS of state depending on system size. 
Other studies related to either sidewards flow or elliptic flow produce similar indecisive results for collision 
energies of the order of a few hundreds MeV [40] . In view of these results we have investigated the sensitivity 
of EE to changes from a soft (A'=210 MeV) to a hard (A'=380 MeV) EoS. Results are shown in left panel 
of Fig. 3; the light and dark bands corresponding to soft and hard EoS respectively; the width of the bands 
were obtained by changing the asy-EoS from a super-stiff {x—-2) to a super-soft {x=2) one. Variations in 
the values of EE due to changes of the compressibility modulus K arc seen to exceed those due to changes 
in the asy-EoS scenario. 

An additional indetermination of theoretical estimates is brought in by the value of the nucleon wave- 
packet width L. Its value is usually set in literature [41] to 2i^=4 fm^ for light systems. Eor heavy systems it 
is found necessary to increase the value to 2L^=8 fm^ in order to generate nuclei with stable static properties 
{e.g. rms). In the right panel of Figure 3 we display results for proton EE for a few values of nucleon wave 
packet width ranging from 2L'^=4 fm^ to 2L^=17 fm^. Variations of the EE magnitude when changing the 
value of L to the one customarily employed for light nuclei or alternatively to the next larger value are 
sizable, of the same order of magnitude as for changes in K from a soft to a hard EoS. 

To summarize, in this section we have shown that while the proton and neutron EEs show important 
sensitivity to changes in the asy-EoS, their equally large of even larger sensitivity to model parameters like 
compressibility modulus K or nucleon wave packet width L limits their usefulness in an attempt to constrain 
the supra-normal density dependence of SE in a model independent fashion. 

3.2. Neutron- Proton Elliptic Flow Difference 

This section is dedicated to an analysis, similar to the one performed for the case of proton EF in the 
previous section, of the neutron-proton elliptic flow difference (npEED): V2~^ = V2 — Due to these 
similarities many of the details of the last Section are identical and therefore have been omitted where 
possible. The interest in this observable stems from the following reasons. The net isospin effect to the total 
EoS is rather small and additionally the pp and nn elementary cross-sections are taken to be the same. One 
would therefore expect that the sensitivity to model parameters of the nucleon elliptic flow is similar 
to the one of protons. An analysis similar to the one performed for protons confirms that to be true. In 
addition, the dependence of the and elliptic flows on the variable x parametrizing the asy-EoS is 
opposite, in line with expectations, due the fact that the neutron asy-potential is repulsive while the proton 
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Figure 3: Left Panel: Sensitivity of the observable to the isospin independent EoS and various choices of asy-EoS. Results 
obtained using a soft (K=210 MeV) or a hard (K=380 MeV) are depicted by light and dark bands respectively. Right Panel: 
Variation of to changes in the width of the wave packet are shown for various values of L (see text). 

asy-potential is attractive. On the experimental side, an accurate measurement of this observable would 
prove challenging as neutron distributions arc rather difficult to measure accurately. 

In Figure (4) a similar analysis as the one performed for protons in Figure 1 but for npEFD is presented. 
The left side panel presents the dependence of this observable on the choice of the asy-EoS. The splitting 
between the extreme 2 and .t = 2, is constant as a function of the impact parameter and its 

magnitude represents about 30-40% of that of either neutron or proton EF the same energy and impact 
parameter. Additionally, the sensitivity to in-medium NN nucleon-nuclcon cross-sections is almost an order 
of magnitude smaller (see right panel of Figure (4)). 

A study of the dependence on model parameters like compressibility modulus and width of the Gaussian 
wave packet shows a sizable suppression as compared to the proton or neutron EF case, see Fig. 5. The 
width of the bands in the left panel of Fig. 5 represents the sensitivity of npEFD to changes of the isospin 
independent EoS from a soft to hard one. The results for various choices of variable x are clearly disentangled 
from each another. Complementary, variations of npEFD due to changes in L, the nucleon wave packet width, 
are also small, of the order of the statistical error of the presented results. 

The results presented in Fig. 4) and 5 for npEFD recommend it as a viable observable for constraining the 
high density behavior of the asy-EoS in a model independent way. The results presented in these figures were 
obtained by simulating Au-f-Au collisions at an incident energy of 400 AMeV and applying the FOPI filter. 
In Figure (6) we extend that study by varying the impact energy from 150 AMeV to 800 AMeV together with 
a variable impact parameter. We plot the impact parameter dependence of A{npE F D)=npF,FD[x = — 2]- 
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Figure 4: Sensitivity of npEFD to the asy-EoS {left panel) and NN in-medium cross-sections {right panel). Statistical errors 
of theoretical estimates amount to about 0.005 in absolute value and have been omitted for clarity. Labeling of curves is the 
same as in Fig. 1. 



npEFD[.T = 2] for various colliding energies, a close to zero value for this variable signaling no sensitivity of 
npEFD to asy-EoS. The most favorable energy domain stretches from T=250 AMeV to almost T=600 AMeV 
and as energy increases in this interval the impact parameter window shrinks to mid-central collisions. For 
incident energies larger than 600 AMeV only mid-central collisions show important sensitivity of npEFD 
to the asy-EoS scenario but with decreasing magnitude as energy increases. The magnitude of the of the 
observed effect demands, for a successful constraining of the asy-EoS, experimental data for the elliptic flow 
of both charged and neutral particles of high quality, of an accuracy of 1% or better. 

We have saved for last a comparison of our estimates for npEFD with the experimental values of neutron- 
hydrogen EFD (nHEFD) [21] and neutron- 'proton' EFD (n'p'EFD) [22] obtained by the FOPI-LAND collab- 
oration. In addition to charged particles, the FOPI-LAND experiment has succeeded to measure azimuthal 
distributions of neutrons, providing experimental evidence of their elliptic flow. Results for charged par- 
ticles with Z—l (hydrogen) have also been published in Rcf. [21]. Extracting pure proton distributions 
is not straightforward as the calorimeter resolution of the FOPI-LAND experimental setup is insufficient 
to separate the H isotopes. In Ref. [22] conservative cuts, that leave out the tail contributions of proton 
energy spectra while retaining some deuteron events, have been applied to extract proton spectra (denoted 
here as 'proton' or 'p' in order to stress the difference from a pure proton spectrum). The FOPI-LAND 
experimental values for proton elliptic flow should therefore be considered as rough estimates only, and due 
to the controversial method employed to extract them a comparison with theoretical estimates should be 
viewed with care. 
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Figure 5: The same as Fig. 3 but for the neutron-proton elliptic flow difference. 



It is well known that flow of heavier isobars is stronger and as a result the values of V2 represent upper 
values (in absolute magnitude) for Wj, nHEFD being therefore an upper value for npEFD. The comparison 
of our estimates for npEFD with the FOPI-LAND nHEFD and n'p'EFT experimental results is presented in 
Figure (7), providing an upper limit to the softness of the asy-EoS. A comparison with theoretical estimates of 
nHEFD should provide tighter constraints; work in that direction is in progress. The n'p'EFD experimental 
data show a scattered dependence on the impact parameter most likely due to the procedure employed 
to eliminate deuteron and triton events from experimental spectra. They are shown here in an attempt 
to stress that further refinements on the experimental side are necessary for a meaningful comparison of 
this observable to theoretical models. The ASY-EOS collaboration, a follow-up of the LAND, FOPI and 
CHIMERA collaborations, is planning to measure both proton and neutron azimuthal distributions with 
an increased accuracy which will have the potential to constrain of the high density dependence of asy-EoS 
through observables like npEFD. 



4. Conclusions 

We have employed an isospin dependent QMD transport model to study the dependence of elliptic 
flow observables on the asy-EoS. The proton and neutron elliptic flows show a monotonous dependence on 
the variable x that parametrizes the various asy-EoS employed, but in opposite directions for protons ans 
neutrons respectively. A sizable dependence of these observable on model parameters like the compressibility 
modulus K, the width of the Gaussian wave packet describing particles in QMD, L, and to a lesser extent 
to in-medium effects on the microscopic NN cross-sections renders them useless for constraining, in a model 

10 



0.1 



I 

II 



> 
I 

II 



> 



0.05 



0.0 



-0.05 




T= 


ISOAMeV 


--- T= 


250 AMeV 


T= 


400 AMeV 


T= 


600 AMeV 


T= 


800 AMeV 



4 

b [fm] 



Figure 6: Energy and impact parameter dependence of the splitting asy supersoft-superstifE of npEFD. 



independent way, the high density dependence of asy-EoS. 

The neutron-proton elhptic flow difference shows (npEFD), in contrast, a sizable dependence on the 
asy-EoS and a almost an order of magnitude smaller dependence on model parameters and NN in-medium 
cross-sections. A study of the dependence of its magnitude on collision energy and impact parameter leads 
to the conclusion that neutron-proton differential squeeze-out can be effectively used to constrain the asy- 
EoS when measured in mid-central heavy-ion collisions at incident beam energies between 250-600 MeV. 
Nevertheless, an effective constraint of the asy-EoS would require experimental values for the elliptic flow 
of both charged and neutral particles of an accuracy of 1% or better. 

A comparison with existing experimental data for the neutron-hydrogen elliptic flow difference, due to the 
FOPl-LAND collaboration, sets and upper limit on the softness of the asy-EoS. The experimental values for 
the neutron-proton elliptic flow chfference show a scattered dependence on the impact parameter, most likely 
due to the controversial method used in the analysis of the experimental data, and as such a comparison with 
theoretical estimates is unreliable and should be viewed with care. The ASY-EOS Collaboration promises 
to remedy these shortcomings and together with a higher accuracy measurement of neutron and proton 
elliptic flows will potentially present the opportunity to constrain the high dependence of asy-EoS through 
observables like neutron-proton elliptic flow difference. 
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